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ABSTRACT 

We present a ^19 ks Chandra ACIS-S observation of the globular cluster Terzan 1. 
Fourteen sources are detected within 1(4 of the cluster center with 2 of these sources 
predicted to be not associated with the cluster (background AGN or foreground ob¬ 
jects). The neutron star X-ray transient, X1732—304, has previously been observed 
in outburst within this globular cluster with the outburst seen to last for at least 12 
years. Here we find 4 sources that are consistent with the ROSAT position for this 
transient, but none of the sources are fully consistent with the position of a radio source 
detected with the VLA that is likely associated with the transient. The most likely 
candidate for the quiescent counterpart of the transient has a relatively soft spectrum 
and an unabsorbed 0.5-10 keV luminosity of 2.6 x 10^^ ergs s“^, quite typical of other 
quiescent neutron stars. Assuming standard core cooling, from the quiescent flux of 
this source we predict long (> 400 yr) quiescent episodes to allow the neutron star to 
cool. Alternatively, enhanced core cooling processes are needed to cool down the core. 
However, if we do not detect the quiescent counterpart of the transient this gives an 
unabsorbed 0.5-10 keV luminosity upper limit of 8 x 10^^ ergs s“^. We also discuss 
other X-ray sources within the globular cluster. From the estimated stellar encounter 
rate of this cluster we find that the n umber of sources we detect is significantly higher 
than expected by the relationship of iPoolev et alJ ((2003)- 

Key words: globular clusters: individual (Terzan 1) — stars: neutron — stars: 
individual (X 1732-304) — X-rays: binaries 


1 INTRODUCTION 

Neutron star X-ray transients form a subgroup of low-mass 
X-ray binaries. Although usually in a quiescent state with 
typical X-ray luminosities of 10®^ - 10®^ ergs s“^, occasion¬ 
ally these systems go into outburst where the luminosity 
increases to around 10®® - 10®® ergs s“®. These outbursts 
are attributed to a large increase in the mass accretion rate 
onto the neutron star. The X-ray spectra from these qui¬ 
escent neutron star transients are usually dominated by a 
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soft component at around 1 keV, and in some cases an addi¬ 
tional power-law component, dominating at energies above 
a few keV, is present. The most widely accepted model used 
to explain the soft, thermal X-ray emission of these tran¬ 
sients in their quiescent states is that in which the emis¬ 
sion is due to the cooling of the neutron star, which has 
been heated during the outbursts. In this case, the quies¬ 
cent lumi nosity should depend on the time-averaged accre¬ 
tion rate llBrown et alJ Il998l : ICamnana et alJ Il998l) . How¬ 
ever, the X-ray emission from the cooling of the neutron 
star cannot directly explain the hard power-law tail and 
its origin is not well understood. Suggestions include resid- 
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ual accretion down to the magnetospheric radius, or pulsar 
shock em ission fe.g..lStella et al.ll994l:ICampana et a,lll 998l: 

ICamnana fc Stellall^nUMenou fc McClintockll2001^. 

Many quiescent neutron stars and other quiescent low- 
mass X-ray binaries have been found in several globular clus¬ 
ters using the Chandra and XMM-Newton X-ray observa¬ 
tions (see iHeinke et al 112003d: iPoolev et a'i]|2003l and ref¬ 
erences therein). Galactic globular clusters provide an ideal 
location to study these types of sources - the distance to host 
clusters can usually be determined more accurately than for 
the Galactic quiescent X-ray binaries. The known distance 
and reddening allows accurate luminosities to be derived and 
removes the distance uncertainty from the quiescent proper¬ 
ties. The high incidence of compact binaries in globular clus¬ 
ters is likely explained by the formation of such binaries via 
exch ange encounters in the very dense environments presen t 
(e.g. IVerbunt fc Hutlll987l : iHut. Murphy fc Verbunti f 199111 . 
though the ultracompact systems may be formed via di rect 
collisions follo wed by orbita l decay lllvanova et al.ll200Rll . 

Terzan 1 llTerzanlll96ll is a globular cluster at a dis¬ 
tance of 5.2 ± 0.5 kpc and a re ddening of E{B — V) = 
2.48 db 0.1 llOrtolani et al.l ^9^. In 1980, the Hakucho 
satellite d etected X-ray bursts f r om a source loca ted in 
Terzan 1 llMakishima et aLlll98ll : llnoue et al.lll98lh . Sev¬ 
eral years later, the persistent source X1732—304 was de- 
tected within Terzan 1 with SL2-XR T onboard Spacelab 2 
JSkhme^^lJ^^S)^ and EXOSAT dWarwick et alJll988t 
IParmar^^nnrSolh This is likely to be the same source as 
the bursting source detected previously by Hakucho. Sub- 

S 3 nt X-ray observation s with ROSAT d.Iohnston et alJ 
; IVerbunt et alJll995^ detected the source with a sim¬ 
ilar luminosity, betwee n 2.0 x 10^^ ergs s"*^ and 1.3 x 10^® 
ergs s“^ (see Fig. 3 of iGuainazzi et S1ll99flll . It is also as¬ 
sumed that X1732—304 is the source of har d X-rays de- 
tected with the SIGMA and ART-P telescopes dBorrel et alJ 
Il996all^ IPavlinskv etlD Il995^ . A possible radio counter¬ 
part of X1732—304 w as observed with the VLA within the 
ROSAT error circles llMarti et al.lll99ljl . 

A BeppoSAX observation of X1732—304 in April 1999 
discovered the source in a particularly low state with the X- 
ray intensity m ore than a factor of 30 0 lower than previous 
measurements llGuainazzi et alJll999^ . A more recent short 
(^ 3.6 ks) Chandra HRC-I observation of Terzan 1 did not 
conclusively detect X1732—304 with a 0.5-10 keV luminosity 
upper limit of (0.5 — 1) x 10^® ergs s~^ depending on the as - 
sumed spectral model dWiinands. Heinke fc Grindlavl2002dl . 
However, an additional X-ray source, CXOGLB J173545.6- 
302900, was detected by this observation. 

In this paper we study the X-ray sources detected in a 
recent ~19 ks Chandra ACIS-S observation of Terzan 1 and 
discuss possible quiescent counterparts of the neutron-star 
X-ray transient X1732—304. 


2 OBSERVATIONS AND ANALYSIS 

Terzan 1 was observed with the Chandra X-ray Observatory 
for ~19 ks on 2005 May 10. The observation was made with 
the Advanced CCD Imaging Spectrometer (ACIS) with the 
telescope aim point on the back-illuminated S3 chip, this in¬ 
creases the sensitivity to low energy X-rays compared to the 
front-illuminated chips. Data reduction was performed us- 



Figure 1. Digitized Sky Survey (DSS) image of Terzan 1. This IR 
plate was originally taken with the UK Schmidt Telescope. The 
inner circle indicates a radius of 0.'7 and the outer circle indicates 
a radius of lf4. 


ing the CIAO v3.2.2 software with the calibration database 
CALDB vS.l.O, provided by the Chandra X-ray Center and 
following the science threads listed on the CIAO website^. 
No background flares were found, so all available data was 
used in our analysis. We removed the pixel randomisation 
added in the standard data processing as this slightly in¬ 
creases the spatial resolution of the images. 

2.1 Source Detection 

The CIAO tool wavdetect ]|2^ was used 

to detect the sources. We detected the sources in the 0.5-8.0 
keV energy band. A detection threshold was chosen to give 
~1 false source over the S3 chip. We detect 39 sources over 
the entire S3 chip with 4 detected counts or greater, after the 
addition of one extra source close to the cluster centre that 
was missed by wavdetect yet clearly a source. Analysing the 
sources that are within the globular cluster half-mass radius 
allows a good balance between including most of the clus¬ 
ter sources and minimising the nu mber of nearby object s or 
backgroun d AGN (see F ig. 1 from IPoolev et al.ll2nn2bL for 
example). iTraeer et alJ lll995^ define the half-mass radius 
for this cluster to be 3.82 arcmin, however, they note that 
their measurement of this value is ‘particularly unreliable’. 
As this value for the half-mass radius is particularly large, 
and given that the majority of X-ray binaries in the cluster 
will be concentrated at the centre we chose a smaller region 
of size 1.'4 within which to analyse the X-ray sources (FiglU 
shows a DSS image of Terzan 1). Within this 1.'4 region we 
find 14 sources. The Chandra image of the cluster is shown 
in Fig. 1^ where all sources within the \'A region are high¬ 
lighted. These sources are listed in descending number of 
detected counts in Table Q also included in this table is a 
list of sources on the S3 chip that are not within \'A of the 
cluster center, listed in decreasing 0.5-8 keV counts. The po¬ 
sition of the brightest source we detect, CXI, is consistent 

^ Available at http://cxc.harvard.edu/ciao/ 
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Figure 2. Chandra ACIS-S image of the globular cluster Terzan 1 in the 0.5-8.0 keV energy band. The large outer circle represents 
the I'.A region within which we analyse the sources. The 14 sources within this region are labelled (in order of decreasing counts in the 
0.5-8.0 keV band) and the extraction regions over laid. Also shown are the la error circles for thr ee ROSAT poin tings (using positions 
for pointings A and B from Ijohnston et ^ 1 ^1995^ and the u pdated position for pointing C from IWrbund 1200511 b the 90% confidence 
level r adius for the EX O S AT observation IParmaret al and the 90% confidence level radius of the radio source detected with the 

VLA iMarti et al.lll99g^ . We predict that 2 of the detected sources are not associated with the cluster. 


with the position of the source CXOGLB J173545.6-302900 
fro m the previous Chan dra HRC-I observation of Terzan 
1 Iw iinands et, 10 Hoi)- From the density of sources on 
the S3 chip outside of the 1.'4 region, we expect to find 2.2 
sources not associated with the cluster (background or fore¬ 
ground objects) within the 1.'4 region (assuming the sources 
outside this region are not associated with the cluster). This 
is in app roximate agreement w ith the log N — log S relation¬ 
ships of ICiacconi et aP (1200 Jl . Thus we expect that 2 of 
the 14 sources we detect within 1^4 are background or fore¬ 
ground sources. 


2.2 Source Extraction 

The IDL tool ACIS Extract iBroos et alJl200^ was used for 
source photometry and extraction of spectra. ACIS Extract 
makes use of CIAO, FTOOLS, ds9 display capability, and 
the TARA IDL software. ACIS Extract constructs polygon 
source extraction regions which are approximate contours 
of the ACIS point-spread function (PSF). The user specifies 
the fraction of the PSF to be enclosed by the contour. The 
shape of the PSF at the location of each source is calculated 
using the CIAO tool mkpsf. For all our sources, except one, 
the PSF fraction contour level was set to 90%, evaluated 
at 1.5 keV. For the brightest source the PSF fraction was 
increased to 95%. With the chosen contour levels there is no 
overlapping of source extraction regions. 
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Table 1. Names, positions (J2000), and background-subtracted counts in four X-ray energy bands are given for the 14 sources detected 
within lf4 of the center of Terzan 1 and the sources detected on the rest of the S3 chip. The statistical error (from wavdetect) in the 
posit ions of all the sourc es is < 0"1, and thus the total positional error is dominated by the uncertainty in the absolute astrometry 
in^b. lAldcroft et al.l2fin7ill . We therefore estimate the Itr positional error for all sources to be ~0”6. The count extraction r egion for each 
source within the 1.'4 region is shown in Fig.|^ The source CXI was first detected (and named) by|Wrinands_et_aJj J200^. 




R.A. 

Dec. 


Net Counts 


Source Name (Label) 

(h m s) 

(O , V ) 

0.5 - 8.0 keV 

0.5 - 1.5 keV 

0.5 - 4.5 keV 

1.5 - 6.0 keV 

(1) 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

CXOGLB J173545.6-302900 (CXI) 

17 

35 

45.57 

-30 29 00.1 

225+1® 

10+3 

157111 

189+1 

CXOGLB J173547.2-302855 (CX2) 

17 

35 

47.26 

-30 28 55.3 

43l® 

1311 

38l^ 

29+ 

CXOGLB J173547.0-302858 (CX3) 

17 

35 

47.02 

-30 28 58.4 

Site 

1511 

3ll6 

16+ 

CXOGLB J173548.3-302920 (CX4) 

17 

35 

48.31 

-30 29 20.8 

29l® 

27+ 

28+ 

+? 

CXOGLB J173545.4-302834 (CX5) 

17 

35 

45.42 

-30 28 34.2 


0+3 

'^-2 

17+ 

15+ 

CXOGLB J173548.7-302940 (CX6) 

17 

35 

48.71 

-30 29 40.0 

nil 

0+2 

11+ 

15+ 

CXOGLB J173549.8-302914 (CX7) 

17 

35 

49.84 

-30 29 14.4 



13+ 

15+ 

CXOGLB J173547.6-302858 (CX8) 

17 

35 

47.67 

-30 28 58.7 

1511 

ll? 

13+ 

14+ 

CXOGLB J173552.4-302931 (CX9) 

17 

35 

52.41 

-30 29 31.3 

7l^ 

0+2 

7+4 

‘ -3 

7+4 

' -3 

CXOGLB J173547.7-302847 (CXIO) 

17 

35 

47.71 

-30 28 47.9 

^-2 


5+3 

^-2 

4+3 

^-2 

CXOGLB J173549.4-302906 (CXll) 

17 

35 

49.42 

-30 29 06.4 

5+3 

^-2 

ll? 

4+3 

^-2 

0+3 

'^-2 

CXOGLB J173543.7-302816 (CX12) 

17 

35 

43.77 

-30 28 16.0 

5+3 

^-2 

ll? 

5+3 

3-2 

4+3 

^-2 

CXOGLB J173547.1-302902 {CX13) 

17 

35 

47.17 

-30 29 02.7 

4+3 

^-2 

0+2 

4+3 

^-2 

4+3 

^-2 

CXOGLB J173542.4-302941 (CX14) 
Sources likely not associated with Terzan 1: 

17 

35 

42.49 

-30 29 41.0 

4+3 

^-2 

ll? 

4 + 

0+3 

'J_2 

CXOU J173538.1-303032 

17 

35 

38.13 

-30 30 32.1 

128+11 

9+3 

^-1 

100 +J 

121+11 

CXOU J173548.1-303420 

17 

35 

48.13 

-30 34 20.2 

7411° 

66l° 

7311° 

8^3 

CXOU J173553.1-303433 

17 

35 

53.12 

-30 34 33.2 

sell 

4611 

571^ 

10+ 

CXOU J173545.6-303204 

17 

35 

45.61 

-30 32 04.9 

4511 

37l^ 

4411 

8l^ 

CXOU J173536.4-302917 

17 

35 

36.49 

-30 29 17.9 

3111 

5+3 

3-2 

30+ 

26+ 

CXOU J173535.9-303143 

17 

35 

35.94 

-30 31 43.4 

3011 

2+? 

21+ 

25+ 

CXOU J173547.1-303426 

17 

35 

47.16 

-30 34 26.1 

2711 

7+4 

'-3 

24+ 

21+ 

CXOU J173557.2-302839 

17 

35 

57.27 

-30 28 39.2 

1911 

0+3 

'^-2 

15+ 

15+ 

CXOU J173559.8-303016 

17 

35 

59.88 

-30 30 16.8 

1711 

5 I 2 

17+ 

12+ 

CXOU J173547.7-303332 

17 

35 

45.78 

-30 33 32.5 

1711 

ll? 

12+ 

14+ 

CXOU J173534.2-302715 

17 

35 

34.28 

-30 27 15.3 

1311 

ll? 

10+ 

11+ 

CXOU J173533.7-303002 

17 

35 

33.79 

-30 30 02.5 

1311 

ll? 

7+4 

' -3 

8+ 

CXOU J173556.0-302717 

17 

35 

56.00 

-30 27 17.3 

1211 

ll? 

7 I 3 

10+ 

CXOU J173550.1-303155 

17 

35 

50.10 

-30 31 55.0 

Oil 

0+3 

'^-2 

10 + 

+3 

CXOU J173554.2-302728 

17 

35 

54.28 

-30 27 28.8 

8ll 

0+2 

8+ 

81^ 

CXOU J173544.2-302736 

17 

35 

44.22 

-30 27 36.2 


0_2 


9+3 

^-1 

CXOU J173603.7-302849 

17 

36 

03.72 

-30 28 49.9 


0+3 

'^-2 


4+3 

^-2 

CXOU J173544.0-303134 

17 

35 

44.07 

-30 31 34.8 


9+3 

^-1 

5+3 

3-2 

5+ 

CXOU J173602.1-303225 

17 

36 

02.11 

-30 32 25.0 

7+4 

‘ -3 

0+2 

7+ 

8+ 

CXOU J173535.9-302650 

17 

35 

35.91 

-30 26 50.9 

7+4 

' -3 

ll? 

6+ 

^-2 

CXOU J173539.7-303007 

17 

35 

39.77 

-30 30 07.0 

^-2 

0+2 

4+3 

^-2 

5+3 

^-2 

CXOU J173555.0-302721 

17 

35 

55.05 

-30 27 21.9 

^-2 

ll? 

5+3 

^-2 

4+3 

^-2 

CXOU J173540.7-302734 

17 

35 

40.72 

-30 27 34.0 

5+3 

^-2 

0+2 

5+3 

^-2 

5+3 

^-2 

CXOU J173547.8-303039 

17 

35 

47.85 

-30 30 39.0 

4+3 

^-2 

+2 

4+3 

^-2 

0+2 

CXOU J173544.7-303125 

17 

35 

44.76 

-30 31 25.6 

4+3 

^-2 

ll? 

4+3 

^-2 

CO Ol 
+ 1 

CO 


For each source a spectrum, lightcurve and event list 
was extracted. The response matrix (RMF) and auxiliary re¬ 
sponse (ARF) files for each source are constructed using the 
CIAO tools mkacisrmf and mkarf. Background subtracted 
photometry was computed for each source in several differ¬ 
ent energy bands: 0.5-1.5 keV (Xsoft), 0.5-4.5 keV (Xmed), 
and 1.5-6.0 keV (Xhard)- These energy bands were chosen to 
be consistent with previous globular c luster studies in the l it- 
erature (e.g. lOrindla.v et al 112001 alibi: IPoolev et al.ll2002allhl: 


iHeinke et alJ I2003allbl3. 120051 : iBassa et alJ l2004^ . Table 0 
lists the counts in each of these bands for all the sources. 
From the counts in the Xgoft, Xmed and Xhard bands we 
create an X-ray color magnitude diagram, plotting the loga¬ 
rithm of the number of counts in the Xmed band against 
the X-ray color, defi ne d as Xcolor = 2.5 l ogXsnft/Xhard 

jGnn^nye^aJ 2001alll3:IPoc^ et al .l2002al lbt lHeinke et alJ 

l2003alrblcl I2005L though we note that Pooley et al. do not 

use the factor of 2.5) (see Fig. |^. Photoelectric absorption 
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Xcolor = 2.5 logio(Xsoft/^hQrd) 

Figure 3. X-ray color magnitude diagram. The lower and 
left axes show the colors and magnitudes from the absorption- 
corrected counts, whilst the upper and right axes show the colors 
and magnitudes determined from the observed counts. We ap¬ 
proximately correct for photoelectric absorption by shifting the 
data -1-0.69 units on the left axis and -^2.25 units on the lower 
axis. Absorption affects each source differently, and our correc¬ 
tion is only an approximation based on typical spectra of globular 
cluster sources. For the purposes of clarity, only a few error bars 
are shown. CX6, CX9 and CX13 are not shown as they have no 
counts in the Xgoft band. 


reduces the number of counts detected and hence affects 
the X-ray colors and magnitudes. To allow comparison with 
sources in other globular clusters this needs to be corrected 
for. However, it is important to note that absorption af¬ 
fects each source differently. We determine an approximate 
correction for the effects of absorption by investigating the 
drop in count rate due to absorption for 3 typical spectra: 
a 3 keV thermal bremsstrahlung, a 0.3 keV bl ackbody, and 
a pow er-law with photon index, F = 2. (e.g. IPoolev et alJ 
l2002all^ . Using PIMMS, we determine the factor by which 
the count rate lowers with the inclusion of photoelectric ab¬ 
sorption at the column density seen f or Terzan 1. The opt i- 
cal reddening, E{B— V) = 2.48±0.1 llOrtolani et al.lll99flll . 
is converted to a column density Nh = 1.36 x 10^^ cm~^ 
(Nh /E(B — V) = 5.5 X 10^^ cm”^, from IPredehl fc Schmitd 
using i? = 3.1). Absorption affects the Xgoft band the 
most, with the absorbed count rate a factor of 13.92 lower 
on average for the three models. The other average correc¬ 
tion factors are 4.85 for X^ed and 1.76 for the Xhard- There 
is not a significant difference in the correction factors for the 
different models. This leads to a shift in Xcolor of 4-2.25 and 
a shift in logXmed of 4-0.69. The lower and left axes of Fig-E] 
show the absorption corrected colors and magnitudes, where 
as the upper and right axes show the observed values. 


2.3 Spectral Analysis 

Using XSPEC (v. 11.3.2) llArnau dll996^ we fit the four bright¬ 
est sources (CXI - CX4) within the U4 region. We attempt 
to fit each of the sources with 3 different models: the neu¬ 
tron star atmosphere model of IZavlin et alJ (Il99(jl . a ther¬ 
mal bremsstrahlung model and a simple power-law. We both 
hx the column density at the cluster value {Nh = 1.36 x 10^^ 


Table 2. Luminosities for X-ray sources in Terzan 1. For CXl- 
CX3 the luminosities have been calculated from the spectral fits 
(see text). For the other sources, PIMMS was used to estimate 
the fluxes from the detected source counts, assuming the cluster 
Nh and a power-law with F = 2.0. No luminosity is given for 
CX4 as it is possibly a foreground object. 


Source 

X-ray Luminosity (ergs 

S 

1 ) 


0.5-2.5 keV 

0.5-10 keV 

CXI 

0.5tJ;3 X 10^3 

2 

^•^-0.3 

X 

1033 

CX2 

1.9tp 7 X 10^2 

^•'^-0.5 

X 

1032 

CX3 

4.3+2;4 X 10^2 


X 

1032 

CX5 

4.9 4i 1.6 X 10^1 


X 

lO^i 

CX6 

X 10^0 


X 

lO^i 

CX7 

2.1+J ® X 1031 

y r^-|-2.3 
'•^-1.8 

X 

1031 

CX8 

l.3tj;7 K 


X 

1031 

CX9 

i-7tj;8 K 

0.i_i 2 

X 

1031 

CXIO 

3.9t3°6° X 1030 

2.6+J® 

X 

1031 

CXll 

S.ltg^s’^ X 1030 

9 9-I-1-6 
^•^-1.2 

X 

1031 

CX12 

l.slo 7 X 1031 

9 9+1-6 

Z.Z_i 2 

X 

1031 

CX13 

3.9tJ°g° X 1030 

y+l-S 

X 

1031 

CX14 

l-7tj;8 X 1031 

y+l-S 

X 

1031 


cm“^) and allow it to be free, using the photoelectric ab¬ 
sorption model phabs. Errors in the fluxes were calculated 
by fixing the free-parameters to their maximum and mini¬ 
mum value in turn, and only one at a time. The model was 
then re-fitted to the data and flux values calculated. Once 
this was done for every free parameter, the flux range was 
used to give the flux errors. 

We discuss each of these sources below. Table El gives 
the luminosities of the Terzan 1 sources as determined by 
spectral fits for CXl-3 assuming a distance to the cluster 
of 5.2 kpc. The luminosities for the other sources are de¬ 
termined by estimating the fluxes from the detected source 
counts using PIMMS, and assuming the cluster Nh and a 
power-law with P = 2.0. A power-law with P = 2.0 cor¬ 
responds to an absorption-corrected Xcolor~ 0.6, which is 
consistent with the majority of these sources (see Fig.EJ. 


2.3.1 CXOGLB J173545.6-302900 (CXI) 

The background-subtracted spectrum was grouped to have 
a minimum of at least 15 counts per bin. Both the neutron 
star atmosphere model and thermal bremsstrahlung model 
do not fit the spectrum satisfactorily with either the cluster 
Nh or allowing it to be a free parameter. A power-law model, 
with absorption at the cluster value, fits with a photon in¬ 
dex, P = 0.2 4z 0.2 and a reduced Cf, = 0.92 for 13 degrees 
of freedom (dof). Allowing the photoelectric absorption to 
be a free parameter leads to a fit with Nh = 2.4)ti ® x 10^^ 
cm“^, P = O.Tj'lQ g, and = 0.82 for 12 dof. The photon 
index of the source from these power-law fits is unusually 
hard for globular cluster sources, and from the X-ray color- 
magnitude diagram it can be seen that this is the hardest 
source in the cluster. Only a few known globular cluster 
sources have P < 1.0 and most tend to have P > 1.5. An 
example of a globular cluster source that under some ob- 
servi ng conditions could look like CXI is the CV X9 in 47 
Tuc llHeinke et al.ll20fCT . This source has very strong soft 
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Figure 4. The spectrum of the brightest source in Terzan 1, 
CXOGLB J173545.6-302900 (CXI). 

X-ray emission (see Fig. 18 in iHeinke et al.l|200Fi) . due to 
oxygen lines. Above ~5 keV it shows more flux than can 
be acconnted for with simple bremsstrahlung or power-law 
models; possibly this is due to partial covering absorption. 
Observing this object only above 2 or 3 keV (for instance, 
if it was in Terzan 1 and heavily obscu red) would give an 
effective photon index of ~0.3. Recently, iMuno et al.l (l2004l 
detected many X-ray sources in the Galactic Center that 
have hard spectra, with a median photon index F = 0.7, 
though many have F < 0. These authors suggest that most 
of the sources are intermediate polars. Intermediate polars 
are more luminous and harder than other CVs. The hard¬ 
ness of these sources is thought to be due to local absorption 
(in addition to Galactic absorption) that is partially cover¬ 
ing the X-ray emitting region, for example by material in 
the accretion flow. We therefore model CXI absorbed by 
the average cluster absorption (using the phabs model, and 
fixing Nh = 1.36 x 10^^ cm“^), plus an additional partially- 
covering absorption component that affects a fraction of the 
emitting region (using the pcfabs model). This leads to a fit 
with xt ~ 0-78 for 11 dof. This spectral fit is shown in Fig.|31 
We get a higher photon index, with F = I.I^q g, with the 
partial absorption component having Nh = 3.4131 x 10^^ 
cm“^ and a covering fraction of Thus, the partial 

covering absorption can explain the extreme hardness as the 
photon index has increased, and the source could be a in¬ 
termediate polar. However, it is not possible to distinguish 
between the quality of these fits given the values of reduced 

xl- 

From this model we calculate the unabsorbed 0.5-2.5 
keV luminosity = O.SlJ'g x 10^® ergs s“^ and the unabsorbed 
0.5-10 keV luminosity = 2.0lJ’3 ^ (assuming 

a distance of 5.2 kpc). 

2.3.2 CXOGLB 3173547.2-302855 (CX2) 

As this source only has a total of 42.8 net counts in the 
0.5-8 keV band, we did not group this spectrum into bins 
and cannot use statistics. We therefore fit the spectrum in 


XSPEC using Cash statistics iCas hll97fl^ without subtracting 
a background spectrum (as Cash statistics cannot be used 
on background subtracted spectra). This should introduce a 
minimal error as there are only 0.2 background photons in 
the extraction region. To investigate the quality of fits, the 
‘goodness’ command is used. This generates 10,000 Monte 
Carlo simulated spectra based on the model. If the model 
provides a good description of the data then approximately 
50% of the simulations should have values of the Cash fit 
statistic that are lower than that of the best-fitting model. 
We define the fit quality as the fraction of the simulated 
spectra with values of the Cash statistic lower than that 
of the data. Both very low and very high values of the fit 
quality indicate that the model is not a good representation 
of the data. 

Fitting a neutron star atmosphere model to this spec¬ 
trum (with the radius = 10 km, the mass = 1.4 Mq and the 
normalisation = l/d(pc)^ = 3.7x10“®) does not give a good 
fit, with a fit quality of 1.0. Fitting a thermal bremsstrahlung 
model also gives a reasonably poor fit, with a fit quality of 
0.85. An improved fit is obtained when this spectrum is fit¬ 
ted with a simple power-law model at the cluster Nh, giving 
a photon index, F = 2.5±0.6 and a fit quality of 0.65. Allow¬ 
ing the Nh to be a free parameter, gives Nh ~ 1 . 2 )')Q'g x 10 ^^ 
cm“^, F = 2.3 ± 1.0 and a fit quality of 0.7. Such a value 
of the photon index suggests that this could be a hybrid 
thermal plus power-law spectrum, but there are not enough 
counts to be conclusive. From the X-ray color-magnitude di¬ 
agram we see that this source is relatively soft and so it is 
possible that it could be a quiescent neutron star X-ray bi¬ 
nary. We determine a 0.5-2.5 keV luminosity of l.Ojjg 7 x 10®^ 
ergs s“^ and a 0.5-10 keV luminosity of 2.6j)Q 5 x 10®^ ergs 
s“^ from the power-law fit with the cluster column density. 

2.3.3 CXOGLB J173547.0-302858 (CX3) 

This source also has too few counts to use x^ statistics and 
so we again use Cash statistics. A simple power-law fit at the 
cluster column density gives a photon index, F = 3.9 ± 1.2 
and a fit quality of 0.51. Such a high photon index strongly 
suggests a soft thermal spectrum, which is also indicated by 
the X-ray color. A neutron star atmosphere fit (with the ra¬ 
dius = 10 km, the mass = 1.4 M© and the normalisation = 
l/d(pc)^ = 3.7x10“®) at the cluster column density gives 

= 73 . 7 )t 3 ® eV with a fit quality of 0.65. Allowing the col¬ 
umn density to be a free parameter does not greatly improve 
the fit, with Nh = l-Slo,® x 10^^ = 78.9 ±6.7 eV 

and a fit quality of 0.62. A thermal bremsstrahlung fit with 
the cluster column density gives kT = 0.6)1)q'| keV with a fit 
quality of 0 . 68 , where as with the column density left free 
we get Nh ~ 1.8!tQ 9 cm“^ and kT = 2 keV with a 

fit quality of 0.77. We determine a 0.5-2.5 keV luminosity 
of 4.312 ® X 10®^ ergs s“® and a 0.5-10 keV luminosity of 
4 . 5 I 2 X 10®^ ergs s“® from the power-law fit with the clus¬ 
ter column density. From its spectrum, this source could be 
a quiescent neutron star. 

2 . 3.4 CXOGLB J173548.3-302920 (CX4) 

Cash statistics have also been used to analyse this source. 
This source is the softest source within the cluster, with an 
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extremely high X-ray color. There is no valid power-law ht 
when hxing the column density to the cluster value. How¬ 
ever, leaving this free we get a fit with Nh = 0.7jlQ g x 10^^ 
cm“^, r = 7 . 1^4 3 and a fit quality of 0 . 74 . There is no 
good fit for either a neutron star atmosphere model or a 
bremsstrahlung model. Such a soft X-ray color and spec¬ 
trum, and the measured Nh being much lower than the 
cluster value suggest that this source may not be associ¬ 
ated with the cluster and maybe a foreground object with 
a lower column density. However, there is no counterpart 
seen in the DSS or 2MASS images. As this could possibly 
be a foreground object and the model is highly uncertain, 
we have not determined a luminosity for this source. 

3 DISCUSSION 

3.1 The neutron-star X-ray transient X1732—304 
in quiescence 

The most likely candidate for the quiescent counterpart of 
the neutron-star X-ray transient X1732—304 is CX2. Of the 

4 sources that are within the ROSAT error circles, it is the 
only one contained in all three and contained within the 
EXOSAT error circle (see Fig. |21 note that the 3 separate 
error circles are from the 3 separate ROSAT observations 
of the source). From the X-ray color-magnitude diagram 
we see that this source is relatively soft and its spectrum 
suggests that it could be a hybrid thermal plus power-law 
source. Thus, it is possible that it could be a quiescent neu¬ 
tron star X-ray binary. The fact that it cannot be fit with a 
po wer-law spectrum alo ne is consistent with what was found 
bv l.Tonker et al.l (|200^; that the power-law fractional con¬ 
tribution increases for sources with luminosities decreasing 
from 10^® erg s“^ to 10®^ erg s“^. The position of CX2 
is consistent wi t h the position of 3 photons detected by 
|w iinands et al.l (l2002^ and has a flux slightly lower than 
their predicted upper-limit. CX3 also has a suitably soft X- 
ray color and spectrum to be a possible candidate, however, 
its position is only fully consistent with the ROSAT A point¬ 
ing and it is on the edge of the error circles for the EXOSAT 
and ROSAT B and C observations. Although CX3 cannot 
be ruled out as the possible quiescent counterpart, CX2 is 
the most likely candidate. 

The source CX2 is also closest to, though not fully con¬ 
sistent with, the position of the radio source detected with 
the VLA which might be the radio counterpart of the tran¬ 
sient. The position of this VLA source (RA: 17^35’”47.'*27, 
Dec: — 30°28'52'.'8) is accurate to 1.7” (90% confidence level 
radius, iMarti et~a DlHi), however, the phase calibrator for 
these VLA observations had calibrator code C, which could 
mean a maximum systematic error of 0715 (J. Marti, private 
communication). The offset between the radio and X-ray po¬ 
sition is 2.5”. The Chandra pointing accuracy is within the 
typica l 0'.'6 absolute astro metry (90% confidence level ra¬ 
dius) iAldcroft et al.ll200(tl as we hnd that the position of 
the X-ray source CXOGLB J173545.6-302900 (CXI) is con¬ 
sistent with the position of the same source in a previous 
Chandra observation of the cluster, with an offset of only 
073. Given the uncertainties in the X-ray and particularly 
the radio position only a small systematic offset is required 
for the sources to be fully consistent. 


For an additional check of the accuracy of the Chan¬ 
dra pointing, possible optical and infra-red counterparts 
were searched for. We searched for USNO-Bl.O (iMonet et alJ 
l2003^ point sources that are close to ACIS Extract posi¬ 
tions for the X-ray sources on the S3 chip. From this we 
get 3 optical sources that match the X-ray positions to 
within the 076 absol ute astrometry (90% conhdence level 
radius) for Chandra JAldcroft et al.l[200(tl . These sources 
are CXOU J173544.2-302736, CXOU J173545.6-303204 and 
CXOU J173536.4-302917 and have offsets of 0707, 0739, and 
0754 respectively. The corresponding USNO-Bl.O sources 
are 0594-0571302, 0595-0586465, and 0594-0571303 which 
have positional accuracies of 072, d!l, 0^'4 respectively (av¬ 
eraging the error in R.A. and Dec.). Matches with the 
2MASS All-Sky Catalog of Point Sources iCutri et alJl200.^ 
were also searched for. Of the USNO-Bl.O matched sources, 
CXOU J173545.6-303204 is the only one to match a 2MASS 
source, with an offset of 0758. CXOU J173550.1-303155 
also has a possible 2MASS counterpart, with an offset of 
Cf.'57. The possible optical (USNO-Bl.O) counterpart for this 
source has an offset of 0778. To check how many accidental 
matches we get, we apply a random offset of up to 20” and 
count the number of matches (sources found to be within 
Cf.'6 of each other). This was repeated 1000 times and the 
number of mean accidental matches of X-ray sources with 
the USNO-Bl.O catalogue was found to be 0.2 ± 0.5. The 3 
sources we hnd is signihcantly above this mean. 

As the radio source is not seen to be fully consistent 
with any X-ray source there is the possibility that it is an 
unrelated source, for example a background AGN or quasar, 
in which case the radio source would likely still be present 
and any future radio observation might prove this. It is also 
possible that we may not have detected X1732—304 in qui¬ 
escence. Only one photon is detected in the 0.5-8 keV band 
within the VLA error circle, which could very likely be a 
background photon. Assuming fewer than 5 counts were de¬ 
tected in VLA error circle, this gives an upper limit on the 
count rate of 2.7 x lO”"* counts s“^. We use XSPEC to esti¬ 
mate an upper limit on the hux from this count rate. Initially 
we Hx the parameters of an absorbed neutron star atmo¬ 
sphere model to the appropriate cluster values expect for 
the effective temperature. This parameter is then adjusted 
to match the required count rate of 2.7 x 10“"^ counts s“^. 
The unabsorbed 0.5-10 keV luminosity upper limit is then 
determined to be 8 x lO®*^ ergs s“^. Such a luminosity is 
sligh tly lower than typically seen for quiescen t neutron stars 
(e.g. iHeinke et al.ll20o'^ : IPooley et al ]|2^ although this 
scenario cannot be ruled out. 

Assuming that CX2 is the quiescent counterpart to 
X1732—304, it is possible to put limits on the length of 
time that the source will stay in quiescence. In the deep 
crustal heating mo del of quiescent neutron star emission 
iBrown et al.ll993) the quiescent luminosity depends on the 
time-ayeraged accretion rate of the source. Assuming stan¬ 
dard core cooling, it is then possible to relate the quiescent 
flux, Fq, and the ayerage flux during outburst, (Fb), with the 
average time the source is in outburst, to, and quiescence, 
fg, v i a Eg » to!{to + tq) X (F’o)/135 fe.g.. IWiinands et alJ 
l20^|2002)- We use the bolometric flux during outburst as 
estimated bv lWiinands et all ll2nn2ll to be 1.5 x 10“® ergs 
cm“^ s“^. To place upper limits on the thermal quiescent 
flux we use the neutron star atmosphere fit to the spectrum 
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of CX2 with the column density at the cluster value. We 
fix the temperature to the upper limit from the fit and de¬ 
termine the bolometric flux (generating a dummy response 
from 0.001 - 100 keV and extrapolating the model). This 
gives Fg < 3.3 x 10“^® ergs cm“^ s“^. The outburst from 
X1732—304 was seen to last at least 12 yr, so, assuming 
that to = 12 yr, we get tg > 395 yr. However, it is not clear 
how long the source was in outburst when it was first de¬ 
tected. If we assume that to — 17 yr, then the lower limit for 
the quiescent time increases to tg > 560 yr. We note that 
the long outburst episode of this transient will have heated 
the crust significantly out of thermal equilibrium with the 
core and therefore the crust might still be cooling down to¬ 
ward equilibrium again. If that is the case, then the thermal 
flux related to the state of the core is even lower than the 
measured value (hence the core is cooler) and the quiescent 
episode should be even longer to allow the core to cool down. 
We note that if CX2 is not the connterpart of the transient 
quiescent flux would be lower and therefore the predicted 
quiescent episodes would be even longer. 

Disk-instability models suggest that such a long quies¬ 
cent should not occur. During quiescence the disk will slowly 
fill up again and at some point, even a small c hange in mas s 
transfer rate should trigger an outburst fe.g.. iLa^ t^ l2nnili . 
Another neutron star X-ray transient observed to have been 
in outburst for >10 yr is KS 1731-260. The quiescent flux 
of this neutron star also leads to the len gth of predicted 

S cent episodes of several hundred years llWiinands et 

. However, the neutron star X-ray transient in the glob¬ 
ular cluster NGC 6440 has much shorter outbursts and qui¬ 
escent pe riods that are consi stent with the observed quies¬ 
cent flux llCackett et aI.ll2'00,'^^ . Thus, it appears that for the 
quasi-persistent neutron stars X1732—304 and KS 1731-260 
to be as cool as observed in quiescence either the quiescent 
periods are several hundreds of years long, or the quiescent 
periods are much shorter and enhanced neutrino cooling is 
required for the neutron star to cool more quickly. In the lat¬ 
ter case it is possible that in these sy stems the neutron stars 
are relatively heavy as suggested bv IColoi et alJ (|20^. 

3.2 Comparison with other globular clusters 

iPoolev et al.l Eiol found a relationship between the num¬ 
ber of X-ray sources (with 0.5-6 keV luminosity Lx > 
4 X 10®° ergs s“®) in a globular cluster and the stellar en¬ 
counter rate of the cluster (see their Figure 2). These authors 
calculate the stellar encounter rate by performing a volume 
integral of the encounter rate per unit volume, R oc p^/v 
(where p is the density and v is the velocity dispersion), 
out to the half-mass radius. This requires knowing the core 
density and velocity dispersion. Unfortunately, for Terzan 1 
there are no measurements of the velocity dispersion avail¬ 
able. However, we can relate the velocity dispersion to the 
core density, po, and core radius, Tc, via the virial theorem, 
and t hus can estimate the encounter rate by JVerbiinfl 

l2nn,3li . As Terzan 1 is a ‘core-collapse’ cluster and the core is 
not well-defined any estimate for the encounter rate based 
on core values will be unreliable. However, using the values 
for the c ore density, p o. and core radius, from the Harris 
catalog iHarri 111993, February 2003 version) we find that 
the core is a factor of 3.4 less dense than that of 47 Tuc and 
the core radius is a factor of 8 times smaller. We therefore 


estimate that the encounter rate of Terzan 1 will be a factor 
of 3.4®’® X 8® ~ 400 smaller than in 47 Tuc. Thus, based 
on the core density and core radius, the encounter rate for 
Terzan 1 is many tens to hundreds of times smaller than the 
value of 47 Tuc. As 47 Tuc has a encounter rate of 396 (in the 
jPoolev et al.l200^ normalisation), we estimat e that Terzan 
1 has an encounter rate of ~ 1, which from the lPoolev et a.lj 
^2003^ relationship would imply that there should only be 
~ 1 X-ray source in Terzan 1. This is clearly not the case, 
as we detect 14 sources (though 2 are predicted to be not 
associated with the cluster). 

As noted above, as Terzan 1 is a core-collapse cluster, 
therefore any estimate based on the core radius and density 
will be unreliable. It is therefore important to determine 
whether any change in these parameters could make this 
cluster consistent with the relationship. If the core radius is 
larger than assumed, this will have the corresponding effect 
of decreasing the computed density of the core, therefore it 
is difficult to increase the estimated encouter rate of Terzan 
1 b y the required facto r of ~50 to make it consistent with 
the lPoolev et al.) (l2003ll relationship. For instance, increas¬ 
ing the assumed core radius by a factor of 3 (to 7'.'2) de- 
crea ses the computed core density by a factor of 3.09 (follow¬ 
ing lDiorgovskilll993|l . and thus will increase the computed 
encounter rate by only a factor of 1.66. Uncertainty in the 
distance has a similarly small effect on the cluster parame¬ 
ters. Uncertainty in the reddening towards Terzan 1 can also 
have an effect on the estimated core density as the core den¬ 
sity is calculated from the extinction correct ed central sur¬ 
face brightness. In fact, using E(B-V)=2.48 iOrtolani et alJ 
I 1999 II . rather than the value from the Harris catalog, will 
increase po by 77%, alth ough still not enoug h to eliminate 
the discrepancy with the lPoolev et al.l (l2003il relation. 

To account for the number of detected X-ray sources 
it may be that the cluster was previously much larger and 
that most of the stars have been lost, perhaps due to pas¬ 
sages th rough the Ga l actic disk as was suggested for NGC 
6397 bv IPoolev et al.l This idea is supported by the 

fact that Terzan 1 is in the bulge of the Galaxy and has the 
closest projection to the Galactic centre among known glob¬ 
ular clusters. Further evidence for a previously more massive 
cluster comes from the spectroscopic study of Terzan 1 by 
lldiart et alJ (l2002fl . These authors identify apparent mem¬ 
bers of Terzan 1 with substantially higher metallicity and 
suggest that Terzan 1 may have captured these stars from 
the bulge during a previous epoch when the cluster was sig¬ 
nificantly more massive. 


4 CONCLUSIONS 

We have presented a ~19 ks Chandra ACIS-S observation 
of the globular cluster Terzan 1. Within 1.4 arcmin of the 
cluster centre we detect 14 X-ray sources and predict that 
2 of these are not associated with the cluster (background 
AGN or foreground objects). The brightest of these sources, 
CXI, is consistent with the position of CXOGLB J173545.6- 
302900 first observ ed during a short Cha ndra HRC-I obser¬ 
vation of Terzan 1 llWiinands et al 11^ - This source has a 
particulary hard spectrum for a globular cluster source, with 
a simple power-law fit giving a photon index of 0.2T0.2. Such 
a hard spectrum suggests that this could be an intermediate 
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polar iIMutio et 3,1.11200^ . The position of the second bright¬ 
est source, CXOGLB J173547.2-302855 (CX2), is the only 
source to have a position that is consistent with all 3 of the 
previous ROSAT poin tings that observed t he neutron star 
transient X1732—304 i.lohnston et a,1.l|l995h as well as the 
EXOSAT observation. This source’s X-ray color and spec¬ 
trum suggest that it could be a quiescent neutron star and 
therefore we find that this source is the most likely candi¬ 
date for the quiescent counterpart of X1732—304. CX2 has 
a position closest to the positi on of the radio so urce de¬ 
tected in Terzan 1 by the VLA llMarti et ahlllOO^i . though 
the positions are not fully consistent. Assuming standard 
core cooling, from the quiescent flux of CX2 we predict ex¬ 
tremely long (>400 yr) quiescent episodes of X1732—304, 
or enhanced core cooling, to allow the neutron star to cool. 
Having estimated the stellar encounter rate of this cluster 
we find sig nificantly mor e sourc es than expected by the rela¬ 
tionship of IPoolev et al.l ^2^^3^ perhaps because the cluster 
was previously much larger and that most of the stars have 
been lost due to passages through the Galactic disk. 
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